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We study the stability of homogeneous magnetic phases in a generalized t — J model includ- 
ing a same-sublattice hopping t' and nearest-neighbor repulsion V by means of the slave fermion- 
Schwinger boson representation of spin operators. At mean-field order we find, in agreement with 
other authors, that the inclusion of further-neighbor hopping and Coulomb repulsion makes the 
compressibility positive, thereby stabilizing at this level the spiral and Neel orders against phase 
separation. However, the consideration of Gaussian fluctuation of order parameters around these 
, , mean-field solutions produces unstable modes in the dynamical matrix for all relevant parameter 

r-{ . values, leaving only reduced stability regions for the Neel phase. We have computed the one-loop 

h-^ ' corrections to the energy in these regions, and have also briefly considered the effects of the correlated 

, hopping term that is obtained in the reduction from the Ifubbard to the t ~ J model. 

'< Two-dimensional antiferromagnets doped with mobile holes have been extensively studied in the last decade in 
^ I order to understand the physics of superconducting cuprates. Most often these studies have been carried out by 
$H . considering the so called t — J model, 



X3 



o 
o 



Hu = ~tJ24a^J. + ^J2^^■S,, (1) 



where i,j are nearest-neighbor sites on a square lattice, the operator cj^ creates a particle on a vacuum where all 
, sites carry a Zhang- Rice singlet, and double occupancy is forbidden (c|^nj — 0). When this Hamiltonian is derived 
Q ' from a realistic three-band Hubbard model for the Cu02 planes, the hopping integral t ~ 2.5 in units of the coupling 
O I J = 1. At half filling (|l|) reduces to the Heisenberg antiferromagnet, which displays quantum Neel order at T = 0. The 
introduction of mobile holes disrupts the magnetic order, and eventually melts it. The description of this quantum 
[ phase transition has lead to the consideration of two related questions: i) the existence of an intermediate spiral 
^ . magnetic phase for small doping, and ii) the possibility of phase separation in the ground state of the t — J model . The 
'■ existence of a spijpl phase is very appealing in view of the incommensurate peaks in the magnetic susceptibility found 
, in real materialsp this incommensurability seems to appear only in superconducting samples, which lends support to 
the idea that is related to hole mobility. Phase separation is possible since for J « t the t — J model is related to 
the large repulsion limit of the Hubbard model, which only accounts for local on-site interactions. Consideration of 
' the true long-range Coulomb interaction, even if strongly screened, could then prevent phase separation. 
] Shraiman and SiggiaH were the first to discuss the existence of a spiral phase in the t — J model. In a semiclassical 
"Y*^ ■ approach, they showed that for infinitesimal doping 5 — 1 — n, where n is the particle density, the lowest-energy spin 
' configuration is a spiral with pitch Q « 5t/ J . Since then, many authors have considered this problem using a variety 
^ ,; of methods. One of the most elegant approaches is the slave fermion-Schwinger boson (SF-SB) theory, in which the 
^ ^ ■ physical particle is replaced by a product of a spinless fermion / and a bosonic spinor b'l' = (6| , foj ) : 

§■ c..^W/, cl^hhl. (2) 

• This replacement, together with the site-occupation constraint J2a4a^io- + fl fi = 1, is a faithful representation of 

1 1 4. 

• 1— I , the original Fermi algebra. Thus, the spins can be written in terms of the Schwinger bosons as Si = ^bl^aa-a-'bia' 
' (cr is the vector of Pauli matrices), while the slave fermion fi accounts for the charge degrees of freedom . This 
^ , representation naturally enforces the t — J model restriction to a single particle per site. Moreover, at half filling 
■ - - ' it correctly reproduaes the spin-wave results for the magnetization and spectrum, unlike the alternative fermionic 
description of spins.B Using the representation (||), Jayaprakash, Krishnamurthy and SarkerQ found a T = phase 
diagram with four phases: a metallic phase with long-range ferromagnetic order for t6 » J, an insulating phase with 
antiferromagnetic correlations for t « J and 5 not too large, a disordered "spin-liquid " insulating phase with only 
short-range antiferromagnetic correlations for t << J and, (5 large, and a spiral metallic phase with incommensurate 



spin order when t5 ~ J. About the same time, Kane et alu argued that the incoherent processes that renormalize the 
hole spectrum in an antiferromagnetic background may be treated by including an effective same-sublattice hopping 
term t' . They considered several homogeneous spin, configurations (Neel, canted, {Q,Q) and {Q,Tr) spirals, double 
spiral) within the same SF-SB mean-field theory oiu. For small doping {5 « 0) these authors found that the spiral 
and Neel orders are favored, with their relative stability depending on the value of t'. They also found that for both 
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5 and t/ J large a state with uniform chiral order is favored. This state is similar to the flux phase obtained near half 
filling in the alternative slave-boson approach, and violates parity and time-reversal symmetries. 

Antiferromagnetic systems are highly susceptible to local deformations of the spin order in the presence of holes. 
Their homogeneous magnetic states are thus inherently unstable due to the large compressibility of the holes, since the 
system wishes to enhance the spin deviations by locally incpeasing the hole density. This tendency to phase separation 
in the t ~ J model was observed by Auerbach and LarsonQ (AL) by studyirig the SF-SB mean-field states by means 
of a expansion in 5. At the mean-field level, the most important result ofOis the negativity of the compressibility k 
for all values of t/J and very small S. This negative k suggests that phase separation between a hole-rich and a Neel 
phase takes place, provided the energy is jininimized by homogeneous states at finite doping. Alternatively|p±he system 
could form local defects or spin polarons.tl Other works using different techniques lead to similar results.l3ii3 However, 
Mori ami Hamadalij showed that the inclusion of the t' term can stabilize a Neel phase for small S. Furthermore, 
Sarkerllil found that in the presence of a Coulomb interaction (screened to a simple nearest-neighbor repulsion V) 
the incommensurate spiral states and also the antiferromagnetic order are favored over the flux states om, and, more 
importantly, flhat they are stable against hole condensation for S larger than some Sc- 

Except fom, all the works mentioned above using the SF-SB representation are mean-field in character. The subtle 
differences among them correspond only to the several ways that the exact form of the constraint is used to simplify 
terms in the Hamiltonian. Since the constraint is later relaxed, all this equivalent Hamiltonians lead to slightly 
different results and, as a consequence, the relative stability of the magnetic phases changes. The more interesting 
question of the general stability of the mean-field picture due to order-parameter fluctuations was partially considered 
by AL, who found that for infinitesimal 5 the spiral phases are unstable for all values of t/J. This instability appears 
as a negative eigenvalue of the dynamical matrix — the Hessian of the ground-state energy as a function of the order 
parameters — . They also considered the effects of the intrasublattice correlated hopping term of 0{t^ /U) that comes 
from the large C//t-jbransformation of the Hubbard model — which was taking into account by the t' hopping term of 
lower symmetry ind — , and stressed the importance of this term in allowing the instability of the magnetic states at 
infinitesimal doping. 

Despite the above reviewed widespread use in the literature of the SF-SB approach, we are not aware of a thorough 
investigation of fluctuations effects on mean-field results. In particular, exactly at half filling {5 — 0) the general SF- 
SB theory reduces to the Schwinger-boson large- iV expansion for the Heisenbpijg antiferromagnet, recently extended 
by usli3 to include the 0{1/N) corrections. In view of the results of Sarker,Eil we believe it is worth reconsidering 
the stability of homogeneous magnetic phases in the t— J — V model since, as stated above, the conclusions of AL 
are restricted to infinitesimal 5 and = 0. In the following we will investigate this question by the consideration 
of order-parameter fluctuations. In particular, the novel features introduced in our calculations are the following: i) 
since the theory presents a local U{1) gauge sjaametry we will use collective coordinate methods — as developed in 
the context of relativistic lattice gauge theoriealj — to handle the infinitely-many zero modes associated to the local 
symmetry breaking in the saddle-point expansion. Contrary to what happens in the radial gauge, in this case the 
dynamical matrix explicitly displays the zero modes, which allows to check the calculations at intermediate steps; ii) 
unlike inEl, here we will consider finite densities and will evaluate the fluctuations corrections in the regionjvhere the 
magnetic phases are stable; iii) we will include a nearest-neighbor repulsion V as considered by SarkerJ12l and also 
a second neighbor hopping t' in the spirit offl. We stress that the.. p i jes ence of a same-sublattice hopping has been 
frequently advocated as a stabilizing mechanism for spiral phases.™ iri The correlated hopping term considered by 
AL will be initially disregarded; its effects will be briefly discussed at the end of this work. We emphasize that these 
calculations are, to the best of our knowledge, the first complete evaluation of fluctuations corrections in the SF-SB 
theory reported in the literature. 

After the Hamiltonian is written in terms of the new operators (|^), the partition function can be expressed as a 
functional integral over coherent Bose and Fermi states. Then, the different terms in ([^) are decoupled by means of 
Hubbard-Stratonovich transformations according to the following scheme (indicated by a vertical bar): 

Ht ^ -Y^t^j hf] \ B.j (3) 
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In Hj we wrote the Schwinger-boson spin-spin interaction in terms of the two SU{2) invariants Bjj = bj^bja- and 

Aij — (jbi^bjff. These operators describe the magnetic fluctuations in the ferromagnetic and antifcrromagnctic 
channels respectively, and allow for a proper treatment of the incommensurate order upon doping. For the hopping 
term Ht we just decoupled bosons from fermions (here the indices i,j run also between next-nearest neighbors). 
Finally, for Hy we indicated the Hartree factorization but we have also considered i) the alternative Fock decoupling, 
and ii) the contribution of this term when the nearest-neighbor repulsion is written in terms of bosons, i.e., Hy = 
yE^J bI I b]^ ; in all these cases the results are qualitatively the same. According to the decompositions ^ of the 
interaction terms, in a saddle-point evaluation of the partition function the Hubbard-Stratonovich fields become the 
following self-consistent order parameters: 1) antifcrromagnctic order: Ai^j ~ (X^o- "'^io-^jo-), 2) ferromagnetic order: 

Bi-j — (Sct^L^jct)) 3) hopping: Fi^j = {flfj), and 4) hole density Ni — (f-fi). The saddle-point Hamiltonian, 
quadratic in boson and fermion operators, can be diagonalized in the standard way. Here we only point out that in 
transforming to momentum space the boson operators we have defined bl^ = J^k b^^^e^'' '^^ with k^- = k + (tQ/2. 
The extra phase aQ.r/2 is required to. obtain the gapless (Goldstone) modes associated to the long-range magnetic 
order at the proper k = 0, ±Q points JIj The self-consistent order parameters are determined using a mesh of 100x100 
points for the integrals in reciprocal space. Since our mean-field theory is not variational because it unphysically 
enlarges the configuration space, we cannot rely only on the energy-minimization criterion to seek for solutions of 
these consistency equations. Consequently, based on previous studies of the t — J model, we considered solutions 
corresponding to magnetic wavevectors (tt — Q, tt) and (tt — Q, tt — Q), which in large regions of parameter space are 
very close in energy. 

At the mean field level we computed the hole compressibility k ^ to investigate indications of possible phase 
separation. In Fig. 1 we plot the energy of the (Q, Q) spiral phase as a function of the density of holes for the standard 
t — J model (t' = = V). From this figure we see that, as found in several studies, the compressibility of holes is 
negative for small densities and all values of t/J, though it becomes positive for large to moderate doping and large 
t/J. The inclusion of a second-neighbor hopping t' makes the spirals and Neel phases to be even closer in energy (less 
than 1% of difference) for all relevant parameter values. In Fig. 2 we give the energy of the Neel phase as a function 
of doping for different values oft' /J and V — 0. We found that this energy does not depend on t, since at mean-field 
order in the Neel phase the dominant processes correspond to coherent hopping in the same sublattice. Fig. 2 shows 
that for t' /J > 0.05 the compressibility becomes positive (the same is valid for the spiral phases and holds in general 
for \t'\/J > 0.05): the non-frustrating character of the second-neighbor hopping promotes hole mobility and prevents 
their segregation in real space. We have also verified that the inclusion of the-Jiearest-neighbor repulsion V further 
favors the stability of the homogeneous magnetic phases, as found by Sarker.O Then, possible instabilities of these 
phases beyond mean-field order are not necessarily related to phase separation. 

We turn now to the consideration of Gaussian fluctuation corrections to mean-field results, which can be obtained 
by extending the calculations described ir£3. The saddle-point expansion breaks the gauge invariance of the theory, 
leading to the existence of saddle-point solutions connected by the continuous U{1) gauge group. This makes the 
dynamical matrix V to have infinitely-many zero modes, which are the Goldstone bosons associated to the (spurious) 
local symmetry breaking. In particular, for translational-invariant saddle-point values of the decoupling fields, there 
is a zero mode 0o(k, uj) — 50{9) /(5^ki^ Ig^g in every k— subspace. Here '^{9) is the vector of gauge-transformed fields, 
which depends on the set of local phases 6, and the (/3^'s can be shown to be right eigenvectors of the nonhermitian 
matrix P. To avoid the infinities associated to these modes without restoring forces we introduce collective coordinates 
along the gauge orbit. Exact integration of these coordinates eliipinates the zero modes and restore the gauge 
symmetry (in the sense that noninvariant operators average to zero) .113 This program can be carried out by enforcittg 
in the functional integral for the partition function the so-called background gauge condition or "natural" gauge,llj 
0{9y .0^ = 0. This condition can be introduced into the functional measure by means of the Faddeev-Popoff trick, 
and restricts the integration to fields fluctuations which are orthogonal to the collective coordinates. At T = 0, after 
carrying out these integrations on the genuine fluctuations, we obtain the one-loop correction to the ground-state 
energy per site, 

27ri_oo V Vl^lv/detI?i(k,L^) / 

Here V± is the projection of T) in the subspace orthogonal to the collective coordinates, and AFp(k, w) = 
\0Q io) .0Q u})\ is the Fadeev-Popoff determinant (i^q (k, w) is the left zero mode of T> in the k — w subspace). 

More explicitly, AFp(k,tj) = [^Y,^ J(l + cosk.r)^^ - 4X^(JS2 + At^B^F^){\ - cosk.r) + uj^] . 

For spiral phases has always negative eigenvalues, which make these phases unstable for all parameter values 
even in the presence of a large repulsion V . Only the Neel phase becomes stable in some regions of the {5, t' / J) plane 
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because of the intrasublattice hopping term, as shown in Fig. 3 for different values of t/ J and V = UnHke the 
mean- field results, in this case the energy depends out/ J since the fluctuations corrections incorporate the incoherent 
intersublattice hopping processes. Fig. 4 shows the corresponding stability diagram for t/ J = 2.5, which is the value 
supposed to be relevant for the cuprates. The consideration of both positive and negative values of t' is interesting 
in view of the proposal of Tohyama and Maekawallll to explain the asymmetry in the behavior of superconducting 
cuprates doped with holes or electrons. According to this proposal, the asymmetry is related to the fact that in the 
generalized t — J model for the Cu02 planes one should have i' > for systems doped with electrons and t' < 
for those doped with holes. The results in Fig. 4 are qualitatively supporting this idea, since the magnetic order 
in systems with t' > (corresponding to electron doping) are more stable than those with t' < (hole doping), as 
observed experimentally. Notice, however, that there is also a region of stability for small negative t' and medium 
to large S, but this could be an artifact of the approximation. This picture is not modified by the presence of the 
nearest-neighbor repulsion V . For completeness, in Fig. 5 we plot the mean-field and fluctuation-corrected energies 
in the Neel phase as a function of doping on the stable part of the line t' = J. 

We have identified the unstable modes in the dynamical matrix. These modes have ~ and are distributed 
over most of the Brillouin zone, although they are more densely located near k = and {tt,tt). They couple the 
fluctuations of the intersublattice fermionic hopping and ferromagnetic order parameter i?k- This clearly shows 
that the tendency of the kinetic term Ht in (^ to locallv, align the spins is responsible for the instabilities of the 
homogeneous magnetic phases, as previously found by AL.eI This effect of local distortions is very difficult to account 
for in mean-field theories, and requires the calculations of fluctuations to treat it. The inclusion of the V term reduces 
the number of unstable modes near k = (as expected, since this term acts against phase separation) but it is not 
able to stabilize the homogeneous phases, a result that had also been conjectured inO. 

Finally, we comment here on the effects of including the three-sitej-|Correlated hopping term of 0{P /U) that can 
be obtained in the reduction of the Hubbard model to the t — J one.B After the replacement (||) this term becomes 
a product of eight boson and fermion operators, requiring a double Hubbard- Stratonovich transformations to be 
decoupled. Using the same order parameters as in (^, we found that its inclusion does not modify qualitatively the 
above results. Ip-particular, it does not play any special roll in stabilizing or destabilizing the magnetic phases. In 
a previous workll3 we considered the alternative Bogoliubov factorization of the four slave-fermion operators in this 
term, which at mean-field order lead to a superconducting phase with some interesting properties, similar to those 
observed experimentally. Unfortunatelly, we found here that when the order-parameter fluctuations are considered 
this superconducting phase becomes unstable like the spiral phases in the magnetic state. 

In this work we have considered, within the SF-SB theory, the general stability of homogeneous magnetic phases in 
a generalized t — J model including a same-sublattice hopping t' and nearest-neighbor repulsion V. At the mean-field 
level we found, in accordance with previous works, that the new terms help to stabilize the Neel and spiral phases 
against phase separation, since for small values of these parameters the compressibility k becomes positive (see Fig. 
2). The correlated hopping term considered in£l does not substantially modify these results. Beyond mean-field theory, 
the study of the dynamical matrix associated to the order-parameter fluctuations reveals that the homogeneous spiral 
phases are unstable in all regions of interest. On the other hand, the Neel phase becomes stable for reasonable values 
of t' > or unphysically large t' < 0, as shown in Figs. 3 and 4. The consideration of moderate values of the 
nearest-neighbor repulsion V does not modify these figures. We stress that Fig. 4 supports qualitatively the proposal 
inim that links the asymmetric behavior of hole-doped and electron-doped cuprates to the t' sign in the corresponding 
effective t — J modeL-for the Cu02 planes. Finally, we found that the consideration of Gaussian corrections to the 
mean-field picture ollla unfortunately destroys the stability of the superconducting phase there obtained. 
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Fig. 1: Mean-field energy E(Q,Q) of the {Q,Q) spiral phase as a function of doping S for different values of t/J 
and t' = = V. 

Fig. 2: Mean-field energy E{7r, tt) of the Neel phase as a function of doping S for different values of t'/J and V = 0. 

Fig. 3: Stability of the Neel phase for t' > 0, V = 0, and different values of t/J. The antiferromagnetic order is 
stable in the regions above the lines. 

Fig. 4: Stability diagram of the Neel phase in the whole {6,t'/J) plane for the physically relevant value t/J = 2.5. 

Fig. 5: Mean- field and fluctuation-corrected energies in the Neel phase as a function of doping S for the parameter 
values t/J = 2.5, t' = J, and V = 0, J. 
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FIG. 4. 




FIG. 5. 
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